arXiv: 1506.04003v2 [cond-mat.str-el] 31Jul2015 


Criticality at the Haldane-insulator charge-density-wave quantum phase transition 


Florian Lange, Satoshi Ejima, and Holger Feliske 
Institut fur Physik, Ernst-Moritz-Arndt-Universitat Greifswald, 17489 Greifswald, Germany 

(Dated: August 4, 2015) 

Exploiting the entanglement concept within a matrix-product-state based infinite density-matrix 
renormalization group approach, we show that the spin-density-wave and bond-order-wave ground 
states of the one-dimensional half-filled extended Hubbard model give way to a symmetry-protected 
topological Haldane state in case an additional alternating ferromagnetic spin interaction is added. 
In the Haldane insulator the lowest entanglement level features a characteristic twofold degeneracy. 
Increasing the ratio between nearest-neighbor and local Coulomb interaction V/U, the enhancement 
of the entanglement entropy, the variation of the charge, spin and neutral gaps, and the dynamical 
spin and density response signal a quantum phase transition to a charge-ordered state. Below a 
critical point, which belongs to the universality class of the tricritical Ising model with central charge 
7/10, the model is critical with c = 1/2 along the transition line. Above this point, the transition 
between the Haldane insulator and charge-density-wave phases becomes first order. 


Topological phases of matter have become one of the 
most fascinating objects of investigation in solid state 
physics HH31- Topological states may arise from topo¬ 
logical band structures or interactions. The order asso¬ 
ciated with these phases can be described by topological 
invariants that are insensitive to gradual changes of the 
system parameters. As a generic feature, topological or¬ 
dered states contain gapless edge excitations that encode 
all the information of bulk topological order [4] . 

Symmetry-protected topological (SPT) phases are 
zero-temperature quantum states with a given symmetry 
and a finite energy gap. The Landau symmetry break¬ 
ing states belong to this class. However, there are more 
interesting SPT states that do not break any symme¬ 
try. For example, in higher dimensions, the Kane-Mele 
band insulator min is a topological state protected by 
U( 1) and time-reversal symmetries. In one dimension, 
a prominent representative is the Haldane phase in the 
spin-1 Heisenberg chain |7], which is protected by inver¬ 
sion, time-reversal, and dihedral symmetries [8j [9]. If at 
least one of these symmetries is not explicitly broken, 
the odd-,!? Haldane insulator (HI) is separated from the 
topologically trivial state by a quantum phase transition. 
Relating topological order and entanglement allows for 
a further classification of topological states [10]. While 
gapped quantum systems without any symmetry split in 
short- and long-range entangled states, the SPT phases 
are always short-range entangled. 

Exploring the connection between topological band 
structures and interacting topological states, it has been 
demonstrated that the topological invariants of gapped 
fermionic systems described by the one-dimensional half- 
filled Peierls-Hubbard model, deep in the Mott insu¬ 
lating regime, can be efficiently computed numerically 
by adding a ferromagnetic spin exchange m On ac¬ 
count of a topological invariant of 2, the Peierls-Hubbard 
model—in a certain parameter regime—possesses the 
same boundary states as the spin-1 Heisenberg chain. 
This raises the question whether the spin-density-wave 


(SDW) and bond-order-wave (BOW) ground states of 
the half-filled extended Hubbard model (EHM) iI2] also 
disappear in favor of a SPT HI phase when a ferromag¬ 
netic spin interaction is added. If the answer is yes, one 
should expect a novel quantum phase transition from the 
SPT state to the charge-density-wave (CDW) insulator. 

In this Rapid Communications, we therefore investi¬ 
gate the ground-state, spectral and dynamical proper¬ 
ties of the EHM with additional, alternating, ferromag¬ 
netic spin coupling J, using the unbiased matrix-product- 
state (MPS) based infinite density-matrix renormaliza¬ 
tion group (DMRG) technique [THHE]- 

The Hamiltonian of the one-dimensional EHM is 

HeHM = ~t^2(4cr C j+la + H - C 0 + U X U o7 n oi 
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Here, cj CT (c Jcr ) creates (annihilates) an electron with spin 
a at site j, rij a = Cj CT c Jcr , t is the transfer amplitude of 
the particles, and U ( V ) denotes their intrasite (intersite) 
Coulomb repulsion. We focus on the half-filled band case. 

The ground-state phase diagram of the EHM has been 
worked out by various analytical na mi and numeri¬ 
cal [HUMS] techniques. In the absence of V (Hub¬ 
bard model), the ground state is a quantum critical 
spin-density wave (SDW) with gapless spin and gapped 
charge excitations V U > 0 [22j . If 2V/U < 1, the 
ground state resembles that at V = 0. For 2 V/U > 1, 
the system becomes a 2 /cf-CDW state, where both spin 
and charge excitation spectra are gapful. The SDW 
and CDW phases are separated by a narrow interven¬ 
ing BOW phase [23, [23] below the critical end point 
[U e ,V e ] = [9.25t,4.76f] [12] . where the ground state ex¬ 
hibits a staggered modulation of the kinetic energy den¬ 
sity (cf. the schematic representations included in Fig. [2] 
below). 
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Here we consider the extended Hamilton operator 

L/2 

H = Heum + J J2 s 2j-iS 2 j (2) 

i=i 

with Sj = (1/2)E<7<t' c L The nearer- 

neighbor Heisenberg spin interaction is assumed to be 
alternating and ferromagnetic, i.e., J < 0 on every other 
bond. Since the EHM at large enough U/V can be 
thought of as spin-1/2 chain, the second term in ([2]) tends 
to form a spin-1 moment out of two spins on sites 2 j — 1 
and 2 j in this limit. Then, the resulting spin-1 antiferro¬ 
magnetic chain may realize a gapped Haldane phase with 
zero-energy edge excitations [9]. 

To proceed we perform an entanglement analysis of 
the model ([2]). The concept of entanglement is inher¬ 
ent in the MPS-based DMRG algorithms too. The so- 
called entanglement spectrum e a characterizes topologi¬ 
cal phases [25], which can be obtained from the singular 
value decomposition. Dividing a system into two sub¬ 
blocks, H = Hl <8> "Hr, and considering the reduced den¬ 
sity matrix p l = Ttr[/j], e a = —21nA a are given by the 
singular values A Q of the reduced density matrix p l. The 
e Q spectrum also provides valuable information about the 
criticality of the system. Adding up the singular values 
A a , we have direct access to the entanglement entropy 
Se = For a critical system with central 

charge c, the entanglement entropy Se between the two 
halves of the infinite chain scales as 

Se = ^ln£x + s o, (3) 

where sq is a non-universal constant. The correlation 
length £ x is determined from the second largest eigen¬ 
value of the transfer matrix for some bond dimension \ 
used in the iDMRG simulation [HHE]. At the critical 
point the physical correlation length diverges, while £ x 
stays finite due to the finite-entanglement cut-off. Never¬ 
theless, £ x can be used to determine the phase transition 
because it increases rapidly with x near the critical point. 
Here, we perform iDMRG runs with x U P to 400, so that 
the effective correlation length at criticality is £ x < 400. 

Let us first discuss the entanglement properties of the 
model ([2]). Figure [l] shows £ x and e a in dependence 
on V/t and U/t for fixed J/t = —1.5. In the weak-to- 
intermediate interaction regime, U/t = 4, we find a pro¬ 
nounced peak in the correlation length at V c /t ~ 2.321, 
which shoots up as x grows from 100 to 200, indicating a 
divergency as £ x —> 00 . Obviously the system passes 
a continuous quantum phase transition. By contrast, 
in the strong interaction regime, U/t = 10, the peak 
height stays almost constant when x is raised. Decreas¬ 
ing |J|, the transition points will approach those of the 
pure EHM, e.g., for J/t = —0.5 we find V c /t ~ 2.242, 
with a simultaneous reduction of the £ x ’s peak heights. 



FIG. 1. (Color online) Correlation length (upper panels) 
and entanglement spectrum (lower panels) as a function 
otV/t for U/t = 4 (left panels) and U/t = 10 (right panels), 
where J/t = —1.5. Data obtained by iDMRG. Dashed lines 
give the BOW-CDW (SDW-CDW) transition for U/t = 4 
(U/t = 10) in the EHM [12]. 

The corresponding entanglement spectra denote that the 
nontrivial phase realized for V < V c resembles the SPT 
Haldane phase of the spin-1 XXZ model [25], in that the 
lowest entanglement level exhibits a characteristic dou¬ 
ble degeneracy [29]. For V > V c , in the CDW phase, this 
level is non-degenerate. 

According to Fig. [T] the maximum in the correlation 
length £ x can be used to pinpoint the HI-CDW quan¬ 
tum phase transition, and with it map out the complete 
ground-state phase diagram of the EHM with ferromag- 



U/t 


FIG. 2. (Color online) iDMRG ground-state phase diagram 
of the one-dimensional (half-filled) extended Hubbard model 
with ferromagnetic spin interaction. The red solid (dotted- 
dashed) lines give the HI-CDW phase boundaries for J/t = 
— 1.5 (—0.5). The quantum phase transition is continuous 
(first order) below (above) the tricritical Ising point [Ut,Vt] 
marked by the star symbol. For comparison the results for 
the BOW-CDW (blue dashed line), SDW-BOW (green dot¬ 
ted line), and SDW-CDW (green double-dotted dashed line) 
transitions of the pure EHM were included m 
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FIG. 3. (Color online) Charge (A c ), spin (A s ), and neutral 
(A n ) gaps as functions of V/t for U/t = 4 (a) and U/t = 10 

(b) . The HI (CDW) phase is marked in grey (white). Panel 

(c) gives the scaling of the entanglement entropy Se(x) with 
the correlation length at the SPT-CDW transition V c /t ~ 
2.321 for U/t = 4. The solid line is a linear fit of the data 
to Eq. <§> indicating an Ising phase transition with c = 1/2. 
Results shown are obtained for J/t = —1.5. 



FIG. 4. (Color online) Central charge c*(L) along the HI- 
CDW transition line for J/t = —1.5. DMRG data (obtained 
with periodic boundary conditions) indicate the Ising univer¬ 
sality class (c = 1/2) for U <Ut and, most notably, a tricrit- 
ical Ising point with c = 7/10 at Ut (red dotted line). Inset: 
Magnitude of the jump of the spin gap as U further increases 
for U >Ut- The infinite MPS data—for a system with infinite 
boundary conditions—point to a first order transition. 


netic spin coupling ([2]). The outcome is given in Fig. [2] 
which also includes the result for the pure EHM (blue and 
green lines). The first striking result is that the HI phase 
completely replaces the SDW and BOW states. That is, 
the HI even survives in the weak-coupling regime untill 
U/t = 0 for any finite J < 0 [provided that V <V C (U, J)]. 
According to this the itinerant model (|2j) behaves as a 
spin-1 model, even at very small U/t where double oc¬ 
cupancy is not largely suppressed. In the intermediate- 
to-strong coupling regime, the HI-CDW transition ap¬ 
proaches the BOW/SDW-CDW transition of the EHM. 
The transition is continuous up to a tricritical Ising point 
[Ut, Vt ](J), which converges to the tricritical point of the 
EHM as J —i 0. In the strong-coupling regime above 
\Ut, Vt], the HI-CDW transition becomes first order. For 
very large U/t the phase boundaries of the HI/SDW- 
CDW transitions are indistinguishable. 

We now characterize the different states and the HI- 
CDW quantum phase transition in more detail. For 
this we first consider the various excitation gaps: A c = 
[E 0 ( N -+ 2,0) + E 0 (N - 2, 0) - 2E 0 ( N, 0)] /2 [(two-particle) 
charge gap], A s = E 0 (N, 1) - E 0 (N, 0) [spin gap], 
and A n = Ei(N,0) — E 0 (N, 0) [neutral gap], where 
E 0 (N e , Sjfot) is the ground-state energy of the finite sys¬ 
tem with L sites for a given number of electrons N e and 
z component of total spin 5f ot , and Ei(N e , Si ot) is the 
corresponding energy of the first excited state. For the 
pure EHM, A c and A n vanish at the BOW-CDW tran¬ 
sition, whereas A s stays finite. Here the excitation gaps 
were determined using DMRG in combination with the 
infinite MPS representation with ‘infinite boundary con¬ 
ditions’ [30U32] . where both finite-size and boundary ef¬ 
fects are significantly reduced. Thereby the whole lattice 
is divided into three parts: a window part, containing 
L\ y sites, and two semi-infinite chains. While the Lw- 


dependence persists, the Lw finite-size scaling is more 
easy to handle than the finite-size scaling in the tradi¬ 
tional DMRG method. Figure [3] shows the variation of 
the different excitation gaps across the HI-CDW transi¬ 
tion in the weak-coupling [Fig. [3](a)] and strong-coupling 
[Fig. [3j[b)] regime. In the former case, the charge and 
spin gaps feature weak minima at the transition point, 
but stay finite. The neutral gap, on the other hand, 
closes, see Fig. [3][a). This is evocative of the Ising transi¬ 
tion between the Haldane and antiferromagnetic phases 
in the spin-1 XXZ model with single-ion anisotropy B8] , 
For U/t = 4, we find V c /t ~ 2.321. In the latter case, 
also the neutral gap stays finite passing the phase tran¬ 
sition [see Fig. [3](b)]. However, the jump of the spin gap 
S s = A S (V/ + ) — A s (V/“) is striking, indicating a first-order 
transition. We obtain V c /t — 5.155 for U/t = 10. 

Next we ascertain the universality class of the HI- 
CDW quantum phase transition. When the system be¬ 
comes critical, the central charge c can easily be deter¬ 
mined from the DMRG entanglement entropy. Utiliz¬ 
ing Eq. ©,Fig . |3j[c) demonstrates that c* indeed follows 
a linear fit to the DMRG data at the critical point (for 
140 < X < 400), provided that prior to that the transi¬ 
tion point was determined with extremely high precision. 
At U/t = 4 and J/t = —1.5, we have c* ~ 0.499(1), 
suggesting the system to be in the Ising universality 
class where c = 1/2. For U/t = 4 and J/t = —0.5 
(not shown), we get c* ~ 0.496(3). From conformal 
field theory [26] the von Neumann entropy for a sys¬ 
tem with periodic boundary conditions takes the form 
Sl(£) = (c/3) ln{(L/7r) sin[(7r£/L)]} + Si with another 
non-universal constant si. With a view of the doubled 
unit cell of the HI phase we slightly modify the related 
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formula for c* |33] : 

* ( 3[S l (L/2-2)-S l (L/2)\ 

1 ln{cos[7r/(L/2)]} ' 1 1 

Figure [4] displays c*{L) when moving along the HI-CDW 
transition line by varying U and V simultaneously. Re¬ 
markably, when U is raised, we find evidence for a 
crossover from c*(L) ~ 1/2 to c*(L) ~ 7/10, which can 
be taken as a sign for an emergent supersymmetry at the 
boundary of the SPT HI phase [M1 - FT(T| . 

Finally, we analyze the dynamical charge (spin) struc¬ 
ture factor of the model ([2]), 

S< z ’\k,w) = Y, l<V>n| Wo}| 2 S(0J - 0J n ) , (5) 


where 6 k = n k (O k = S z ). In Eq. (J5), \^ 0 ) (| 4> n )) de¬ 
notes the ground (nth excited) state, and ui n = E n — Eq. 
Following Ref. [JO], we first evaluate the related two- 
point correlation functions, ( , ipo\Oj(T)Oo(0)\'ipo), by way 
of real-time evolution of the ground-state infinite MPS 
IV’o) ■ Thereby we apply infinite boundary conditions to 
a finite window of sites (L\y = 128). After O is applied to 
a given site the system is evolved at least up to r = 30/t, 
where a time step St = 0.05/f is used in the fourth-order 
Suzuki-Trotter decomposition. Fourier transformation 
then gives the dynamical structure factors. 

For the spin-1 chain and extended Bose-Hubbard mod¬ 
els it has been demonstrated that the dynamical spin and 
density structure factor reveal distinguishing features in 
the SPT and topologically trivial phases [23 [37J S8]. Fig¬ 
ure [5] illustrates the intensity of the dynamical wave- 
vector-resolved spin and density response in the k-u 
plane. In the HI phase, both S zz {k,uj) and S(k,uj) ex¬ 
hibit an essentially symmetric line shape with respect to 
k = 7t/2 and gaps at k = 0 and 7r, but the spectral weight 
of the excitations is higher for k > 7r/2; see Figs.[5])a) and 
§b). While the spin response remains unaffected at the 
Ising transition point [Fig. [5](c)] , the gaps in the charge 
response closes at k = n, reflecting the doubling of the 
lattice period CDW phase [Fig. [5jd)] . Obviously S(k,ui) 
follows the behavior of the neutral gap rather than those 
of the charge gap [cf. Fig. [3](b)] . In the CDW phase, the 
overall lineshape of S zz {k : hj) is asymmetric with a larger 
excitation gap at k = n. Note that we find now two 
dispersive features (branches) in S zz (k,u) and S(k,Lo), 
where a changeover of the intensity maximum takes place 
at k = 7t/2 [cf. Fig. [3](e) and (f)]. 

To summarize, exploiting the link between topologi¬ 
cal order and entanglement properties, we examined the 
ground-state and spectral properties of the paradigmatic 
one-dimensional extended Hubbard model (EHM) with 
alternating ferromagnetic spin coupling J by numerically 
exact (DMRG) techniques. We showed that any finite 
spin interaction J < 0 stabilizes a symmetry-protected 
topological Haldane insulator (SPT HI) that replaces the 
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FIG. 5. (Color online) Intensity plots of the dynamical spin 
structure factor S zz (k,ui) (top) and density structure factor 
S{k,u) (bottom) in the SPT HI phase (left), at the HI-CDW 
transition point (middle), and in the CDW phase (right). 
Dashed lines connect the intensity maxima at given fc. Other 
model parameters are U/t = 4 and J/t = —1.5. 


spin-density-wave and bond-order-wave ground states ex¬ 
isting in the pure EHM below a critical ratio of nearest- 
neighbor ( V ) to intrasite (U) Coulomb interaction. The 
HI manifests the twofold degeneracy of the lowest entan¬ 
glement level and, regarding the dynamical spin/density 
response, reveals a similar behavior as the SPT state of 
the spin-1 chain [28] and the HI of the extended Bose- 
Hubbard model Emus. Furthermore, analyzing the cor¬ 
relation length, entanglement spectrum and many-body 
excitation gaps, we found clear evidence for a quantum 
phase transition from the SPT HI phase to a CDW when 
the V/U -ratio is raised. Using iDMRG, the HI-CDW 
boundary and therefore the complete ground-state phase 
diagram could be determined with very high accuracy. 
In the weak-to-intermediate interaction regime, the HI- 
CDW transition belongs to the Ising universality class. 
Here the central charge c = 1/2, and only the neutral gap 
vanishes. This is reflected in the dynamical density struc¬ 
ture factor, where the gap closes at momentum k = n, 
just as for the Hl-antiferromagnet transition of the spin-1 
chain. In the strong interaction regime we found a first- 
order phase transition characterized by a jump in the spin 
gap. Decreasing the magnitude of J, the HI-CDW phase 
boundary approaches the BOW-CDW transition line in 
the pure EHM; thus, making the system topological, this 
changeover can be determined more precisely. Perhaps 
most interesting, tracing the central charge along the HI- 
CDW transition line, we detect a tricritical Ising point 
with c = 7/10 that separates the continuous and first- 
order transition regimes. A further field theoretical study 
would be highly desirable to elucidate the origin of the 
tricritical Ising point. In either case the EHM with ad- 
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ditional ferromagnetic spin exchange provides valuable 
insights into the criticality and nontrivial topological ex¬ 
citations of low-dimensional correlated electron systems. 
Note that we applied the ferromagnetic spin exchange in 
order to easily realize an effective spin-1 state. Including 
a physically more relevant dimerization of the transfer 
intergrals (hopping) will also stabilize the HI phase, so 
that the Ising quantum phase transition occurs between 
the HI and CDW phases [JS). Then, in this extended 
Peierls-Hubbard model, the tricritical Ising point with 
c = 7/10 will separate the HI-CDW transition line into 
continuous and first-order lines m- 

Note added in proof. Due to the quantum-classical cor¬ 
respondence D-dimensional quantum and (D + l) dimen¬ 
sional classical systems share important physical proper¬ 
ties. So it is well known that the quantum spin-1 chain is 
related to the classical two-dimensional restricted-solid- 
on solid (RSOS) model [3T]. It has been shown that a 
Fibonacci anyonic chain can be mapped-using the RSOS 
representation of the algebra-onto the tricritical Ising 
model with c = 7/10. The transitions observed in our 
model can be understood as transitions from a low- 
density phase to a high-density phase of doubly occu¬ 
pied sites. Interestingly in the hard squares model first 
and second order transitions from low to high densities 
also occur with a tricritical point with the same central 
charge c = 7/10, see Ref. [42] and references therein. This 
connects-at the tricritical point-our model, hard squares 
and the so-called golden chain [35] , 
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